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Abstract. We address the dilepton production from the parton interactions in the early 
stage of relativistic heavy-ion collisions within the parton-hadron-string dynamics (PHSD) off- 
shell transport approach. The description of partons in PHSD is based on the dynamical 
quasiparticle model (DQPM) matched to reproduce lattice QCD results in thermodynamic 
equilibrium. According to the DQPM the constituents of the strongly interacting quark-gluon 
plasma (sQGP) are massive and off-shell quasi-particles (quarks and gluons) with broad spectral 
functions. In order to address the electromagnetic radiation of the sQGP, we derive off-shell 
cross sections of — >■ 7* , — > 7* -f g and qg j*q {qg 7*g) reactions taking into account 
the effective propagators for quarks and gluons from the DQPM. Dilepton production in In+In 
collisions at 158 AGeV and in Au-f Au at ^/s = 200 GeV is calculated by implementing these 
off-shell processes into the PHSD transport approach. By comparing to the data of the NA60 
and PHENIX Collaborations, we study the relative importance of different dilepton production 
mechanisms and point out the regions in phase space where partonic channels are dominant. 



Dileptons are emitted over the entire space-time evolution of the heavy-ion colhsion, from 
the initial nucleon-nucleon collisions through the hot and dense phase and to the hadron decays 
after the freeze-out. This is both a challenge and advantage of the probe. The separation of 
different "physics" in the dilepton radiation is nontrivial due to the non-equilibrium nature of the 
heavy-ion reactions, and covariant transport models have to be used to disentangle the various 
sources that contribute to the final dilepton spectra seen experimentally. 

1. The PHSD approach to dilepton radiation 

To address the dilepton production in a hot and dense medium - as created in heavy-ion 
collisions - we employ an up-to-date relativistic transport model, the Parton Hadron String 
Dynamics [T] (PHSD), that incorporates the explicit partonic phase in the early reaction region. 
Within PHSD, one solves generalized transport equations on the basis of the off-shell Kadanoff- 
Baym equations for Green's functions in phase-space representation (in the first order gradient 
expansion beyond the quasiparticle approximation). The approach consistently describes the 
full evolution of a relativistic heavy-ion collision from the initial hard scatterings and string 
formation through the dynamical deconfinement phase transition to the quark-gluon plasma 
(QGP) as well as hadronization and to the subsequent interactions in the hadronic phase. In the 
hadronic sector PHSD is equivalent to the Hadron-String-Dynamics (HSD) transport approach 
[21 [3l S] that has been used for the description of pA and AA collisions from SIS to RHIC 
energies. In particular, (P)HSD incorporates off-shell dynamics for vector mesons - according 



to Refs. [5j - and a set of vector-meson spectral functions [6j that covers possible scenarios for 
their in-medium modification. 

In the off-shell transport description, the hadron spectral functions change dynamically 
during the propagation through the medium and evolve towards the on-shell spectral function 
in the vacuum. As demonstrated in [6] the off-shell dynamics is important for resonances with a 
rather long lifetime in vacuum but strongly decreasing lifetime in the nuclear medium (especially 
oj and (j) mesons) and also proves vital for the correct description of dilepton decays of p mesons 
with masses close to the two pion decay threshold. 

Let us remind the off-shell transport equations (see [7] for details). One starts with a first 
order gradient expansion of the Wigner transformed Kadanoff-Baym equation using the following 
expression for the memory integrals 

/><(xi,x2)= -/^^d^x'G(ti-t') [S>(xi,x')-S<(xi,x')] G><(X',X2) (1) 

+ f^^d^x' T.><{xi,x') Q{t2-t') [G>{X',X2)-G<{X',X2)] 
= -I-ood^'x' S^(xi,x') G><ix',X2) + S><(xi,x') G^(x',X2) . (2) 

and arrives at a generalized transport equation [5l |8] : 

2p''d=' i G>< - {S"^ + Re S^, iG><} -{i S><, Re G^} = i iG> - i i:> i G< (3) 



{M, i G><} -{i S><, Re G^} = i^< iG> -ii:> i G< (4) 

as well generalized mass-shell equation 

[p^ -m'^ -Y.^ -Re ^^]iG>< = iT>< Re G^ + \ {iS>,iG<} - \ {iE<,iG>} (5) 
U 

with the mass-function M. In the transport equation ([3]) one recognizes on the l.h.s. the drift 
term p^^ d^^ iG^^, as well as the Vlasov term with the local self-energy T,^ and the real part of 
the retarded self-energy Re S^. On the other hand the r.h.s. represents the collision term with 
its typical 'gain and loss' structure. 
The mass-function for fermions is 

Mf{p,x) = ul-U'-ml\ (6) 

with the effective mass and four-momentum given by 

^hi^^P) = mh + U^{x,p), U''{x,p) = p'' -Ul^{x,p) , (7) 

where U^{x,p) and Uf^{x,p) denote the real part of the scalar and vector self-energies of the 
particle, respectively, and stands for its the bare (vacuum) mass. After inserting ([6]) into 
the generalized transport equation ([3|), the covariant off-shell transport theory emerges, that has 
been denoted as HSD [H d] . It is formally written as a coupled set of transport equations for 
the phase-space distributions Nh{x,p) (x = {t,f), p = {e,p)) of fermion /i [2ll3j with a spectral 
function Ah{x,p) (using iG^{x,p) = Nh{x,p)Ahix,p)), i.e. [Qj 

{(n^ - n^dPUi: - mldlUi) + {u^d-^Uj; + mld'^U^) d^} Nh{x,p) Ah{x,p) 

-{i^<,Re G«} = Eh.h.h, Tr2Tr3Tu [T^T]i2^uS\U + Ha - - n4) (8) 

xAh{x,p)Ah2{x,P2)Ah,_^{x,P3)Ah^{x,P4,) 
X {Nh^{x,P3)Nh^{x,P4)fh{x,p)fh2{x,P2) - Nh{x,p)Nh2{x,P2)fh3{x,P3)fh4{x,P4)} ■ (9) 
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Figure 1. Diagrams contributing to the dilepton production from a QGP: (a) Born mechanism, 
(b) gluon-Compton scattering (GCS), (c) vertex correction, (d) gluon Bremsstrahlung 
(NLODY), where virtual photons (wavy hues) spht into lepton pairs; spiral lines denote gluons, 
arrows denote quarks. In each diagram the time runs from left to right. 

Here = {dt, Vr) and df, = {de, Vp), (/i = 0, 1, 2, 3). The backflow term in ^ is given by 

- {^^<,ReG^} = ( M,(..p1^+^rf(l,p)V4 ) i^^^^^P) (10) 

-^M {iM^wM^) i^'^i^^p) r/.(^,p)) • (11) 

It stands for the off-shell evolution which vanishes in the on-shell limit, when the spectral 
function Ah{x,p) does not change its shape during the propagation through the medium, i.e. 
VrT{x,p)=0 and VpT{x,p)=0. The set of coupled differential-integral equations ([8]) is solved 
using a test-particle ansatz. 

The transport description of quarks and gluons in PHSD is based on the knowledge of the 
partonic propagators from a dynamical quasiparticle model for partons matched to reproduce 
lattice QCD results in thermodynamic equilibrium (DQPM). The DQPM describes QCD 
properties in terms of single-particle Green's functions (in the sense of a two-particle irreducible 
approach) and leads to the notion of the constituents of the sQGP being effective quasiparticles, 
which are massive and have broad spectral functions (due to large interaction rates). The 
transition from partonic to hadronic degrees of freedom in PHSD is described by covariant 
transition rates for the fusion of quark-antiquark pairs to mesonic resonances or three quarks 
(antiquarks) to baryonic states. 

Dilepton radiation by the constituents of the QGP proceeds via the elementary processes 
illustrated in Fig. [TJ the Born (Drell-Yan) q + q annihilation mechanism, Gluon Compton 
scattering (q + g ^ ^* + q and q + g ^ ^* + q), and quark + anti-quark annihilation with 
gluon Bremsstrahlung in the final state {q + q ^ g + 'J*)- In the on-shell approximation, one 
uses perturbative QCD cross sections for the processes listed above. However, in order to make 
quantitative predictions at experimentally relevant low dilepton mass and strong coupling, we 
have to take into account the non-perturbative spectral functions and self-energies of the quarks, 
anti-quarks and gluons thus going beyond the on-shell approximation. For this purpose, off-shell 
cross sections were derived in [10] for dilepton production in the reactions q + q ^ (Drell- 
Yan mechanism), q + q — )• g + (quark annihilation with the gluon Bremsstrahlung in 
the final state), q{q) + g ^ q{q) + l^l^ (gluon Compton scattering), g ^ q + q + and 
9(9) ~^ q{q) + 5 + (virtual gluon decay, virtual quark decay) in the sQGP in effective 
perturbation theory by dressing the quark and gluon lines with the DQPM propagators for 
quarks and gluons. The obtained off-shell elementary cross sections then are implemented into 
the PHSD transport code, where the masses of quarks and gluons are distributed according to 
the DQPM spectral functions. 

For illustration, we plot the dilepton production cross sections in the Born mechanism in 




Figure 2. Dilepton production cross sections 
in the Born channel {q + q ^ + ^~)- The 
short dashed (black) line shows the on-shell, 
i.e. the standard pQCD result. The other lines 
show the ofF-shell cross section, in which the 
annihilating quark and antiquark have finite 
masses mi and m2 with different values. 



Figure 3. The HSD results for the mass 
differential dilepton spectra in case of inclusive 
Au + Au collisions at ^/s = 200 GeV in 
comparison to the data from PHENIX \18\ 
IT9] in the 'free' scenario, the 'collisional 
broadening' picture as well as the 'dropping 
mass + collisional broadening' model. 



Fig. [2j The short dashes (black) line shows the on-shell, i.e. the standard perturbative 
result. The other lines show the off-shell cross section, in which the annihilating quark and 
antiquark have finite masses mi and m2 with different values: mi = 0.3 GeV, m2 = 0.15 GeV 
(solid magenta line), mi = 0.3 GeV, m2 = 0.3 GeV (dash-dotted red line), mi = 0.3 GeV, 
m2 = 0.6 GeV (dashed blue line), mi = 0.6 GeV, m2 = 0.6 GeV (dash-dot-dot green line). 

2. Comparison to data 

By employing the HSD approach to the low mass dilepton production in relativistic heavy-ion 
collisions, it was shown in [JJJ that the NA60 Collaboration data for the invariant mass spectra 
for f^^n^ pairs from In+In collisions at 158 A-GeV favored the 'melting scenario pL2j. On the 
other hand, the dilepton spectrum from In+In collisions at 158 A-GeV for M > 1 GeV could 
not be accounted for by the known hadronic sources (see Fig. 2 of [llj). Also, hadronic models 
do not reproduce the softening of the mr distribution of dileptons at M > 1 GeV [T2] . 

In Refs. [131 El ES] we presented PHSD results for the dilepton spectrum as produced in 
In + In reactions at 158 AGeV compared to the NA60 data [12\ I16j . We confirmed the HSD 
results that the spectrum at invariant masses below 1 GeV was better reproduced by the p 
meson yield, if a broadening of the meson spectral function in the medium was assumed. On the 
other hand, the spectrum at M > 1 GeV was shown to be dominated by the partonic sources. 
Moreover, accounting for partonic dilepton sources allowed to reproduce in PHSD the effective 
temperature of the dileptons (slope parameters) in the intermediate mass range. The softening 
of the transverse mass spectrum with growing invariant mass implies that the partonic channels 
occur dominantly before the collective radial flow has developed. 

The PHENIX Collaboration has presented dilepton data from pp and Au + Au collisions at 
Relativistic-Heavy- Ion-Collider (RHIC) energies of -^=200 GeV [T71 [18l [19] , which show a large 
enhancement in Au + Au reactions (relative to scaled pp collisions) in the invariant mass regime 
from 0.15 to 0.6 GeV/c^. Moreover, if realistic partial loss of the D- and l)-meson correlations 
due to their rescattering is kept in mind, one has to conclude that there exists another domain 




Figure 4. The PHSD results for the mass 
differential dilepton spectra in case of inclusive 
Au + Au collisions at ^/s = 200 GeV in 
comparison to the data from PHENIX [181 US] 
in the low mass region (M = — 1.2 GeV). 



Figure 5. The PHSD results for the mass 
differential dilepton spectra in case of inclusive 
Au + Au collisions at ^/s = 200 GeV in 
comparison to the data from PHENIX [18\ [T9] 
for M = 0-4 GeV. 



of invariant mass, in which the measured dilepton yield in A + ^ collisions is underestimated by 
the scaled yield from p + p; this domain is at masses from 1 to 4 GeV/c^ \14\ I20j. 

Let us first present in Fig [3] the HSD results for e~^e~ pairs in inclusive Au + Au collisions in 
comparison to the data from PHENIX [18. ,19j as calculated in [llj , recalling that HSD provides 
a reasonable description of hadron production in Au + Au collisions at y/s = 200 GeV [21 J. 
Whereas the total yield is quite well described in the region of the pion Dalitz decay as well 
as around the co, (f) and J/^ mass, HSD clearly underestimates the measured spectra in the 
regime from 0.2 to 0.6 GeV by approximately a factor of 5 [TT]. After including the in-medium 
modification scenarios for the vector mesons, we achieve a sum spectrum which is only slightly 
enhanced compared to the 'free' scenario (see Fig. 5). The low mass dilepton spectra from 
Au + Au collisions at RHIC (from the PHENIX Collaboration) are clearly underestimated in 
the invariant mass range from 0.2 to 0.6 GeV in the 'collisional broadening' scenario as well as 
in the 'dropping mass + collisional broadening' model. We mention that HSD results for the 
low mass dileptons are very close to the calculated spectra from van Hees and Rapp as well as 
Dusling and Zahed [22] (cf. the comparison in Ref. \19\ I23j). At higher masses (from 1 to 4 
GeV) the only hadronic sources of correlated lepton pairs are the charmed mesons: semi-leptonic 
decays of correlated D-mesons and the dilepton decays of charmonia. Between the (p and J/^ 
peaks, the HSD results underestimate the PHENIX data by approximately a factor of two. 

By implementing the off-shell partonic processes into the PHSD transport approach, we 
calculate the dilepton spectra in Au + Au at ^/s=200 GeV and compare to the PHENIX data 
in Figs, m and [5l In Fig. [J] we show our results for low masses (M = — 1.2 GeV); in this 
region, the yield in PHSD is dominated by hadronic sources and essentially coincides with the 
HSD result. Note that the collisional broadening scenario for the modification of the /)- meson 
was used in the calculations presented in the Figs. [l]and[5j There is a discrepancy between the 
PHSD calculations and the data in the region of masses from 0.2 to 0.6 GeV. The discrepancy 
is not amended by accounting for the radiation from the QGP, since the latter is 'over-shone' 
by the radiation from hadrons integrated over the evolution of the collision. In contrast, the 
partonic radiation is considerable in the mass region M = 1 — 4 GeV as seen in Fig. [5j The 
dileptons generated by the quark-antiquark annihilation in the sQGP constitute about half of 



the observed yield in the mass range between the masses of the (f) and the J/^* mesons. For 
M > 2.5 GeV the partonic yield dominates over the D-meson contribution. Thus, accounting 
for partonic radiation in PHSD fills up the gap between the hadronic model results |1H [20 ] and 
the data at Af > 1 GeV. 

3. Summary 

The Parton Hadron String Dynamics [Ij (PHSD) transport approach incorporates the relevant 
off-shell dynamics of the vector mesons as well as the explicit partonic phase in the early hot 
and dense reaction region. A comparison of the transport calculations to the data of the NA60 
Collaborations points towards a 'melting' of the p-meson at high densities, i.e. a broadening of 
the vector meson's spectral function. On the other hand, the spectrum for M > 1 GeV is shown 
to be dominated by the partonic sources. 

The low mass dilepton spectra from Au + Au collisions at RHIC (from the PHENIX 
Collaboration) are clearly underestimated by the hadronic channels in the invariant mass range 
from 0.2 to 0.6 GeV. The discrepancy is not amended by accounting for the radiation from the 
QGP, since the latter is 'over-shone' by the radiation from hadrons integrated over the evolution 
of the collision. In contrast, the partonic radiation is visible in the mass region M = 1 — 4 GeV. 
The dileptons generated by the quark-antiquark annihilation in the sQGP constitute about half 
of the observed yield in the mass range between the masses of the (p and the J/^ mesons. For 
M > 2.5 GeV the partonic yield dominates over the D-meson contribution. Thus, accounting 
for partonic radiation in PHSD fills up the gap between the hadronic model results and the data 
for M > 1 GeV. 
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